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Abstract: Exploring the surface energy exchange between atmosphere and water bodies is essential to 
gain a quantitative understanding of regional climate change, especially for the lakes in the desert. In this 
study, measurements of energy flux and water vapor were performed over a lake in the Badain Jaran 
Desert, China from March 2012 to March 2013. The studied lake had about a 2-month frozen period 
(December and January) and a 10-month open-water period (February-November). Latent heat flux (LE) 
and sensible heat flux (Hs) acquired using the eddy covariance technique were argued by measurements of 
longwave and shortwave radiation. Both fluxes of longwave and shortwave radiation showed seasonal 
dynamics and daily fluctuations during the study period. The reflected solar radiation was much higher in 
winter than in other seasons. LE exhibited diurnal and seasonal variations. On a daily scale, LE was low in 
the morning and peaked in the afternoon. From spring (April) to winter (January), the diurnal amplitude 
of LE decreased slowly. LE was the dominant heat flux throughout the year and consumed most of the 
energy from the lake. Generally speaking, LE was mostly affected by changes in the ambient wind speed, 
while Hs was primarily affected by the product of water-air temperature difference and wind speed. The 
diurnal LE and Hs were negatively correlated in the open-water period. The variations in Hs and LE over 
the lake were differed from those on the nearby land surface. The mean evaporation rate on the lake was 
about 4.0 mm/d over the entire year, and the cumulative annual evaporation rate was 1445 mm/a. The 
cumulative annual evaporation was 10 times larger than the cumulative annual precipitation. Furthermore, 
the average evaporation rates over the frozen period and open-water period were approximately 0.6 and 
5.0 mm/d, respectively. These results can be used to analyze the water balance and quantify the source of 
lake water in the Badain Jaran Desert. 
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1 Introduction 


Lakes cover a large area of the Earth's surface and make vital contributions to regional energy 
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balance and the water cycle. Because the surface energy exchange, with different heat capacity 
and albedo from those of the surrounding land surface, through biogeochemical and biophysical 
processes, lakes can affect climate to some extent at both local and regional scales, even at the 
global scale (Bonan, 1995; Small et al., 1999; Magnuson et al., 2000; Walter et al., 2006; Long et 
al., 2007; Liu et al., 2012; Lee et al., 2014). Environmental changes, open-water surface 
hydrological processes, and water—atmosphere system can all have a significant influence on 
surface energy exchange, and it is important to gain a quantitative understanding of how these 
factors affect climate and weather (Hostetler and Bartlein, 1990; Bates et al., 1993, 1995; Tsuang 
et al., 2001; Li et al., 2016). Some previous studies on lakes have found that lakes can act as 
substantial energy reservoirs that introduce massive seasonal thermal lags in the surrounding 
landscape due to their large heat capacities of lakes (Blanken et al., 2000; Schertzer et al., 2003; 
Downing et al., 2006; Nordbo et al., 2011). Intraseasonal surface evaporation and energy radiation 
of lakes are significant and are closely corrected with strong winds and other weather events 
(Rouse et al., 2003, 2005; Lenters et al., 2005; Liu et al., 2009; Kettle et al., 2012). However, 
there have been few studies focusing on the surface evaporation and energy radiation of lakes in 
regions with strong winds, especially the desert regions. 

Eddy covariance (EC) techniques have been widely and successfully used to measure carbon, 
water and energy fluxes between atmosphere and water bodies, forests and other terrestrial 
ecosystems (Barford et al., 2001; Xu et al., 2013; Baldocchi, 2014; McGloin et al., 2014; Jung et 
al., 2017), as well as lakes (Blanken et al., 2000; Liu et al., 2009; Xiao et al., 2013; Biermann et 
al., 2014). However, few EC measurements have been made over inland lakes; moreover, little or 
no work has been done on lakes in deserts because of the harsh environment. Meanwhile, it 
should be noted that the energy flux depends on many aspects, such as lake size, lake depth and 
surrounding terrain. As a consequence, EC measurements need to be conducted in regions with 
different climate and weather conditions (Lenters et al., 2005; Liu et al., 2009). Furthermore, the 
EC observation network needs to be supplemented by studies of lakes in deserts to provide 
reliable and comprehensive knowledge of energy and water vapor fluxes over different terrains. 

The Badain Jaran Desert is one of the driest regions in China. However, more than 110 
perennial lakes are distributed among megadunes in its hinterland (Wang et al., 2016). These 
lakes vary in shape, size and salinity. Previous studies have investigated the geochemistry of the 
groundwater and lakes (Yang and Williams, 2003; Chen et al., 2004; Gates et al., 2008a, b), the 
lake-level fluctuations in the past recorded by palaeoshorelines (Yang et al., 2010), the aeolian 
geomorphology and evolution of megadunes (Yang and Williams, 2003; Zhu et al., 2012; Rioual 
et al., 2013), and the key scientific problem, namely, the source of lake water (Dong et al., 2016). 
From the perspective of water balance, the water vapor exchange is an important part of water 
budget; and from the perspective of energy balance, the lake-atmosphere exchange is a significant 
part of energy transfer, especially the evaporation of lakes. Accurately quantifying the evaporation 
of lakes is a particular research challenge, especially in the water-scarce area (Xu et al., 2013; 
McGloin et al., 2014). Study of the water cycle and energy flux in the Badain Jaran Desert is of 
great significance for understanding the lake-atmosphere interaction and the source of lake water, 
as well as for managing water resource. 

The aims of this study are quantifying the energy flux over a lake in the Badain Jaran Desert to 
provide detailed knowledge about the thermal conditions, estimating the evaporation of the lake 
using EC measurements, and determining the relationships of sensible heat flux (Hs) and latent 
heat flux (LE) with several environmental factors. This study may fill a gap in the current 
knowledge of lake-atmosphere interactions over lakes in arid regions, especially in deserts. The 
results could provide a scientific reference for future studies on the source of lake water. 


2 Materials and methods 


2.1 Study area 
The Badain Jaran Desert (39°04'15"—42°12'23"N, 99°23'18”—104°34'02"E) lies in the western 
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part of Inner Mongolia Autonomous Region, China. It is surrounded by the ancient Juyan Lake 
and Guaizi Lake to the north, the Beida Mountain and Heli Mountain to the south, the ancient 
Gurinai Lake and the Zhengyi Gorge of the Heihe River to the west, and the Yabrai Mountain to 
the southeast (Ma et al., 2014). The Badain Jaran Desert has an area of about 52x10? km? (Zhu et 
al., 2010). The prevailing wind direction in the interior and surrounding regions of the desert is 
northwesterly in spring, easterly in summer and autumn, and westerly in winter (Chen, 2011). The 
wind speed in the northern part of the desert is much greater than that in the southern part of the 
desert (Yang et al., 2003; Dong et al., 2004). The annual mean wind speed varies from 2.8 to 4.6 
m/s, and the highest wind speed is observed in April and May of each year. The annual mean 
diurnal temperature is 34.4°C (Ma et al., 2014), and the mean temperatures in winter and summer 
are —9.1°C and 25.3°C, respectively. The Badain Jaran Desert is therefore classified as a "cold 
desert” (Warner, 2004). The mean annual precipitation is 90.1-115.4 mm in the south of the desert 
and 35.2-42.9 mm in the north of the desert. Thus, the climate is characterized by a continental 
arid climate. Shallow groundwater and springs are present around lakes in the desert. The total 
dissolved solids (TDS) of the springs are less than 1 g/L. The lakes are primarily located in the 
southeast of the desert, i.e., on leeward slopes of dunes in the interdune depressions. 

Yindeertu Lake, located in the Badain Jaran Desert, is a large lake with a small island far from 
the shore (39°50'N, 102°27'E; 1169 m a.s.l.; Fig. 1). It has a surface area of approximately 1.03 
km?. The measured maximum depth is approximately 9.4 m, and the mean depth is about 5.6 m. 
No runoff has been found over the lake. The TDS content in this lake exhibits significant seasonal 
variations, ranging from 140 to 220 g/L. During the study period (March 2012 to March 2013), 
freeze-up and ice break-up occurred in late December and early March, respectively. The 
dominant plant species around the lake are Artemisia desterorum and Glycyrrhiza yunnanensis. 
The lake is surrounded by dunes, with heights ranging from 200 to 300 m. The dunes and lakes in 
the Badain Jaran Desert are relatively stable and do not move within one year. 


2.2 Experimental measurements 


The eddy covariance (EC) measurement system (E1) was located in the center of the Yindeertu 
Lake, 250 m from the nearest shore (Fig. 1c). The space is enough to satisfy the flux measurement 
for all wind directions. The radiometers were mounted on the water surface beside the EC system. 
The data for this study were collected from March 2012 to March 2013. 

In principle, the EC measurement method provides direct estimates of momentum flux, 
sensible heat flux (Hs) and latent heat flux (LE). In this study, these turbulent fluxes were 
measured by an EC system, which included a three-dimensional sonic anemometer (R3-50, GILL, 
UK) and an open-path CO/H2O infrared gas analyzer (LI-7500A, Li-COR, USA). These 
instruments were set at the height of 2.3 m and were separated by a horizontal distance of 20 cm. 
Wind velocities in three dimensions and atmospheric temperature fluctuations were measured, 
together with variations in water vapor density. Sensor signals were acquired with a logger 
(LI-7550, Li-COR, USA) at 10 Hz. Relative humidity and air temperature were obtained using 
Humicap sensors (HMP155, Vaisala, Finland) installed at the height of 3.0 m on airtight 
double-steel radiation shields. The water surface temperature was determined using an infrared 
radiometer (Apogee, SI-111, USA), which was installed at the height of 1.5 m above the water 
surface beside the EC system. The net radiation was obtained through a four-component net 
radiometer (NRO1, Hukseflux, Netherlands) at the height of 1.5 m above the water surface, as 
with the infrared radiometer. A tipping-bucket rain gauge (HOBO RG3-M, Onset, USA) installed 
at the height of 1.5 m was employed to record the total precipitation at 30-min intervals 
automatically. Slow-response sensor signals were obtained with a CR-3000 data logger 
(Campbell Scientific Inc., USA), at an average interval of 30 min. All the sensors were 
fast-response and small enough to cause no flow disturbances, meeting the requirements for the 
EC technique (Lorrai et al., 2010; Lemaire et al., 2017). 

Another EC measurement system (E2) was located on a nearby flat sandy land on the lower 
part of a megadune beside Sumubarunjilin, 6.68 km away from the EC on the lake (E1). This EC 
system included the same types of instruments as the lake-based EC system (Hu et al., 2015). 
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Fig. 1 Distribution of the eddy covariance (EC) instrument on the Yindeertu Lake (E1) and land surface (E2) in 
the lake group, southeastern Badain Jaran Desert (a), location of the E1 on the Yindeertu Lake (b), and landscape 
of the Yindeertu Lake in the Badain Jaran Desert (c) 


2.3 Data processing 


The monthly averages of air temperature, water surface temperature, albedo over the lake, wind 
speed, relative humidity and water vapor density, and the total precipitation were obtained 
through the EC system. 

The values of Hs (W/m?) and LE (W/m?) were determined using the EC method as follows 
(Eqs. 1 and 2, respectively): 


Hs=p,c,wT , (1) 


LE=p,L,wq, (2) 
where pa (kg/m?) is air density; cp (J/(kg-K)) is the specific heat of air at a constant pressure, i.e., 
1004 J/(kg-K); L, (J/kg) is the latent heat of water vapor, which can be obtained with an equation 
of air temperature (Aubinet et al., 2000); and the parameters w’ (m/s), 7’ (K) and q’ (kg/kg) are the 
deviations from the temporal averages of the vertical wind speed, air temperature and specific 
humidity, respectively. Furthermore, the overbars in Equations 1 and 2 indicate the time average 
(over 30 min in this study). 
The component fluxes of the radiation balance were expressed as Equation 3: 
Rn=DR-—UR-(ULR-DLR), (3) 
where Rn, DR and UR are the net all-wave radiation, incident solar radiation and reflected solar 
radiation, respectively; ULR is the outgoing longwave radiation from the lake surface; and DLR 
is the incoming longwave sky radiation. The units of all these parameters are W/m. 

Monthly variations of Rn, DR, UR, DLR, ULR, ULR-DLR, the daily total evaporation rate 
(ET) and the seasonal and daily variations of LE and Hs were explored, in four months 
representing different seasons: January (winter), April (spring), July (summer) and October 
(autumn). Correlation analyses were performed to obtain the driving forces behind heat fluxes. 
Finally, comparisons of heat fluxes between the lake and land surface, and between evaporation 
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and precipitation were obtained. 
2.4 Data quality control and flux footprint analysis of EC 


After periodic calibration and instrument maintenance, the raw data were obtained from the EC 
system at 10 Hz. These data were calculated using the EddyPro Express post-processing software 
(Version 6.2, Li-COR, USA) in a recommended manner. All data were processed with a rigorous 
quality control through a uniform procedure. 

The key step of the EC technique involves both calculation and correction. After much debate 
and examination of the results from field observations, we established standards for data 
collection, calculation and correction. Specifically, we removed the raw data outliers and applied 
corrections to the raw data for angle-of-attack errors. These errors were caused by imperfections 
in the sonic anemometer (Gash and Dolman, 2003; Van der Molen et al., 2004) and were dealt 
with using Nakai correction functions (Nakai et al., 2006). A double rotation was applied to the 
outliers (Aubinet et al., 2000). A compensation was used to take account of the time lag between 
the gas analyzer and sonic anemometer measurements. The average vertical deflection in the wind 
direction was —1.1°. The frequency response correction (Moncrieff et al., 2005), sonic 
temperature correction (Van Dijk et al., 2004) and density fluctuation correction (WPL-correction) 
(Webb et al., 1980) were applied (Hu et al., 2015). We also carried out a turbulence stationarity 
test (Foken and Wichura, 1996) and an overall turbulence characteristic test, evaluated the quality 
of data and divided the data according to the quality level. Furthermore, we determined data that 
were needed to be interpolated or that were missing. The interpolated data accounted for 29% of 
the total samples. 

There were several examples of gap filling (Aubinet et al., 2000; Falge et al., 2001; Liu et al., 
2013). It should be noted that data recorded during dew events and rainy days were not utilized to 
avoid errors in the calculation resulting from the presence of liquid water on the sensor window of 
the CO2/H20. An important thermodynamic factor, i.e., the energy of evaporation, mainly came 
from the DR, as did Hs and LE. The data on DR were relatively reliable and continuous. 
Therefore, gaps in daytime data were filled using the relationship between the DR and the 
measured Hs, LE and ET. In this study, the data gaps were filled by cumulative calculation. If the 
data were indeed missing for less than 3 d, they were interpolated. In contrast, data missing for 
more than 3 d were not interpolated unless there was a need for accumulated data. Data gaps at 
nights were filled with data from the following and preceding dates by linear interpolation. The 
half-hourly Hs and LE were finally obtained using the EC method. 

It is important to quantify the sampling area of the turbulent flux measurements (Hs and LE), 
which is also known as the flux footprint and is defined as the upwind area that contributes to the 
flux. We assessed the flux footprint using the Lagrangian particle dispersion model proposed by 
Kljun et al. (2004). During the daytime, there was no apparent prevailing wind direction on the 
lake, although the northerly frequency was relatively large. However, there was an obvious 
change in prevailing wind direction between day and night: southerly during the daytime and 
northwesterly during the nighttime. The source areas making the most important contribution 
over the whole day, during the nighttime and during the daytime were respectively at 40.6, 38.2 
and 42.9 m distances from the EC system. About 70% of the source areas were at 137.3, 138.5 
and 136.2 m distances from the EC system over the whole day, during the nighttime and during 
the daytime, respectively. Because the distance from the nearest lake shore was greater than 250 
m, these results indicated that the observation conditions substantially met the demands for 
successful observation of fluxes due to the data were obtained from the lake. 


3 Results 


3.1 General characteristics of the over-lake meteorology 


The monthly averages of air temperature during the observation period showed that from March 
2012 to March 2013, the highest air temperature was found in July 2012, with a monthly average 
temperature of 27.7°C, while the lowest air temperature was found in January 2013, with a 
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monthly average temperature of —6.7°C. The air temperature over the lake was above 0°C from 
March to December and below 0°C from December to February next year (Fig. 2a). The 
maximum (35.5°C) and minimum (—19.3°C) half-hourly averages of temperature were found in 
August and December, respectively (Fig. 2a). Based on the observed surface temperature and 
albedo over the lake from March 2012 to March 2013, we concluded that the lake had about a 
2-month frozen period (December and January). We took these two months (i.e., December and 
January) as a basis for comparison. When the lake was frozen, the variations in lake surface 
temperature were quite similar to those in air temperature. The highest and lowest lake surface 
temperatures were 30.1°C and —13.4°C, respectively (Fig. 2b). The mean wind speed was 1.7 m/s 
at the height of 2.5 m, and the maximum half-hourly average wind speed was 10.1 m/s. The wind 
speed in summer was greater than that in winter (Fig. 2c). The relative humidity was fluctuated 
during the study period, with the average value of 35.9% (Fig. 2d). It should be noted that no 
humidity data were collected in August because of the malfunction of the sensor. The total 
precipitation was 145 mm during the study period, with most of them occurring from March to 
September (Fig. 2e). There was an extreme precipitation event occurred on 20 July 2012 (47 mm 
precipitation), which was extremely heavy in the desert area. 
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Fig. 2 Half-hourly averages of air temperature (a), water surface temperature (b), wind speed (c) and relative 
humidity (d), and cumulative precipitation (e) from March 2012 to March 2013. The data of air temperature, 
water surface temperature, wind speed and relative humidity were smoothed over a 48-h period. It should be 
noted that no humidity data were collected in August because of the malfunction of the sensor. 


3.2 Variations in radiation 


Variations of Rn, DR, UR, DLR, ULR, and ULR—DLR from March 2012 to March 2013 are 
shown in Figure 3. Generally speaking, Rn, DR, DLR and ULR from the lake surface exhibited 
significant seasonal variations (Figs. 3a, b, d and e). However, UR exhibited a completely 
different pattern (Fig. 3e). It was significantly higher in January and February than in other 
months. This may be due to that the amount of UR was closely related to the underlying surface 
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conditions. The lake was covered with ice from December 2012 to February 2013. The albedo of 
the ice surface was much higher than that of the open water. Therefore, the value of UR was high 
during this period. Note that there was no snow on the lake during the frozen period and thus the 
albedo was completely from the ice without any contribution of snow. 
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Fig. 3 Variations of net all-wave radiation (Rn, a), solar radiation (DR, b), reflected solar radiation (UR, c), 
incoming longwave sky radiation (DLR, d), outgoing longwave radiation (ULR, e) and ULR—DLR (f) from 
March 2012 to March 2013 


Seasonal and daily variations in the five types of radiation (including Rn, DR, UR, DLR and 
ULR) are presented in Figure 4. The diurnal amplitude of the longwave radiation (i.e., DLR and 
ULR) was smaller than that of the shortwave radiation (i.e., DR and UR). 

The net longwave radiation on the surface of the lake (i.e., ULR-DLR) was always positive 
during the study period (Fig. 3f), which indicated that the lake continuously lost energy through 
longwave radiation. The net longwave radiation exhibited two peaks in June and November. 

Both shortwave and longwave radiation fluxes showed seasonal dynamics and daily 
fluctuations; however, they differed specifically at several aspects. DR exhibited smooth and 
symmetric changes during the day and peaked at solar noon (Fig. 4b), suggesting that clouds are 
relatively rare in the Badain Jaran Desert. In comparison, the peak value of UR was larger in 
January than in other months, and UR exhibited more complex and non-symmetric changes in all 
months with the exception of January (Fig. 4c). The values of ULR and DLR were positive during 
both the daytime and nighttime, with the daytime values being greater than the nighttime values 
(Figs. 4d and e). The fluctuations of ULR and DLR were smaller than those of the other types of 
radiation (including Rn, DR and UR). The variations of Rn were similar to those of DR, with the 
positive peaks occurring at noon (Figs. 4a and b). Furthermore, the values of Rn were negative at 
night, and the largest amplitude occurred in midsummer. 
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In general, all the radiation values were larger in summer (July) than in winter (January), 
except for the UR values, which were greater in winter than in summer (Fig. 4). The maximum 
Rn and DR values usually occurred during 12:00—14:00 (LST), the maximum UR value appeared 
at 09:30 in October, and the maximum DLR and ULR values were observed between 15:00 and 
18:00. 
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Fig. 4 Seasonal and daily variations of net all-wave radiation (Rn, a), solar radiation (DR, b), reflected solar 
radiation (UR, c), incoming longwave sky radiation (DLR, d) and outgoing longwave radiation (ULR, e) in April, 
July, October and January 


3.3 Variations in heat fluxes 


During the study period, LE exhibited remarkable seasonal changes, whereas Hs did not show 
obvious seasonal changes (Fig. 5). LE increased from January to July and decreased from July to 
January next year. The minimum and maximum values of LE were —25.3 and 799.0 W/m?, 
respectively, with an average value of 122.7 W/m”. The positive LE values mean that heat is lost 
from the lake via evaporation, while the negative LE values suggest the occurrence of 
condensation. Over the course of the year, LE was predominantly positive, and negative LE only 
occurred occasionally in winter. The average Hs was 5.7 W/m’, which was much smaller than the 
average LE. Hs varied slightly throughout the year, with large values being observed in winter. 
Compared to Hs, LE exhibited clear seasonal and diurnal variations (Fig. 6). On a daily scale, 
LE was low in the morning and peaked in the afternoon, whereas the maximum Hs occurred in 
the morning and the minimum Hs occurred in the afternoon or the evening. The results suggested 
that the available energy was primarily consumed by latent heating in the afternoon and by 
sensible heating in the morning. From spring (April) to winter (January), the diurnal amplitude of 
LE decreased slowly, while the changes in the diurnal amplitude of Hs varied slightly. Rn 
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Fig.5 Variations of sensible heat flux (Hs) and latent heat flux (LE) from March 2012 to March 2013 
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Fig. 6 Seasonal and daily variations of sensible heat flux (Hs), latent heat flux (LE) and net all-wave radiation 
(Rn) in April (a), July (b), October (c) and January (d) 


exhibited obvious difference between the daytime and nighttime periods, with the values being 
positive during the daytime and negative during the nighttime. The loss of Rn was mainly via the 
LE throughout the year (except for winter). Rn might be influenced by wind in spring and autumn 
and by solar radiation in summer. 
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As shown in Figure 7a, the daily total ET peaked in summer. From 1 April 2012 to 31 March 
2013, the maximum ET was 11.74 mm/d in June, and the mean ET was about 4.0 mm/d. 
Furthermore, the average ET was approximately 0.6 mm/d during the frozen period (December 
and January) and 5.0 mm/d during the open-water period (February—-November). Approximately 
40% and 60% of evaporation occurred during the nighttime and daytime, respectively. Moreover, 
ET generally decreased throughout late summer and autumn. It maintained relatively low values 
from December to January and increased from March to June. The average ET was lower in July 
than in June and August. It should be noted that the maximum ET did not occur in July, although 
the temperature was highest in this month. The cumulative evaporation was 1445 mm from 1 
April 2012 to 31 March 2013 (Fig. 7b), which was 10 times larger than the accumulated 
precipitation over the same period. 
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Fig. 7 Daily total evaporation rate (ET, a) and cumulative evaporation (b) from 1 April 2012 to 31 March 2013 


4 Discussion 


4.1 Radiation and heat fluxes 


Annual variations in DR were largely affected by solar elevation and weather conditions, 
including reflection, scattering and absorption due to ambient atmospheric conditions. Monthly 
variations in DR were due to the changes in the azimuthal angle of the sun, while daily variations 
were due to the effects of clouds. The UR was determined mainly by the features of the 
underlying surface conditions, particularly the albedo. The studied lake was frozen in winter 
(January), and the albedo was higher in winter than in other seasons. Therefore, UR was relatively 
large in winter (Fig. 4c). The diurnal curve of the albedo was not symmetric. This may be due to 
the highly variable sand hills surrounding the lake, which tended to block incident UR. The 
monthly average albedo was 0.13 from March to November; however, it was 0.21 from December 
to February, which is greater than the value (i.e., 0.07) reported for a water reservoir in a 
semi-arid region (Gallego-Elvira et al., 2010) but smaller than the value in an alpine meadow or a 
savanna (Beringer and Tapper, 2002; Al-Riahi et al., 2003; Berbert and Costa, 2003; Zhang et al., 
2010). However, the albedos obtained in this study are similar to the average values of the Earth 


SUN Jie et al.: Eddy covariance measurements of water vapor and energy flux over a lake in the... 


and the Northern Hemisphere (Gupta et al., 1999). This suggests that the energy loss from the 
lake ecosystem via radiation may be lower than the loss from other ecosystems because of the 
relatively low albedo. 

The exchange of longwave radiation between atmosphere and lake surface is controlled 
primarily by surface features, clouds, atmospheric temperature, and atmospheric water vapor 
pressure. Over the studied lake, the water vapor pressure varied from 0.1 to 2.0 KPa (Fig. 8). 
However, it was previously reported that the water vapor pressure in most ecosystems typically 
varied from 2.0 to 5.0 KPa (Kellner, 2001; Hunt et al., 2002; Wever et al., 2002). The water vapor 
pressure in the studied desert lake ecosystem was slightly higher than that in the deserts (ranging 
from 0.01 to 1.5 KPa) and much smaller than that in most other ecosystems. Because of the low 
ambient water vapor pressure and relatively few clouds, DLR over the studied lake was obviously 
lower than the value given by NASA for the Earth as a whole (Gupta et al., 1999). 
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Fig. 8 Variations in water vapor pressure over the lake (a) and the surrounding land (desert) surface (b) from 
March 2012 to March 2013 


The diurnal amplitude of LE was smaller in winter and larger in spring, which indicated that 
the ice cover effectively restricted the heat exchange between air and water. Simultaneously, the 
ice cover significantly increased the surface albedo and reduced the absorbed incident shortwave 
radiation. 


4.2 Driving forces behind heat fluxes 


Turbulent exchanges of water vapor and heat flux over lakes are affected by the surrounding 
environment. The heat flux over a lake surface is determined mainly by the ambient temperature, 
turbulent mixing intensity, and water-air vapor pressure deficit (Hostetler and Bartlein, 1990; 
Bonan, 1995; Liu et al., 2009; Nordbo et al., 2011). The air at the interface between the lake and 
the overlying atmosphere is always saturated, and the saturation pressure is a function of the lake 
surface temperature at the interface (Hostetler and Bartlein, 1990). Meanwhile, several other 
parameterizations have been proposed, especially for lake evaporation (Rosenberry et al., 2007). 
After calculating the relationship between Hs and LE by taking account of several 
environmental factors (Fig. 9), we draw similar conclusions as provided by previous studies. 
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Specifically, there were obvious correspondences between Hs and water-air temperature 
difference (Nordbo et al., 2011; Liu et al. 2015), between Hs and wind speedxwater-air 
temperature difference (Blanken et al., 2000, 2003; Liu et al., 2009; Nordbo et al., 2011), between 
LE and wind speed (Liu et al., 2009), and between LE and wind speedxwater-air vapor pressure 
deficit (Blanken et al., 2003; Liu et al., 2009; Nordbo et al., 2011). According to the study of 
Blanken et al. (2000), LE could be better described by wind speedxwater-air vapor pressure 
deficit than by water-air vapor pressure deficit alone at Great Slave Lake. Taking account of the 
influence of wind on latent heat and evaporation, we analyzed the correlations of LE with wind 
speed and wind speedxwater-air vapor pressure deficit. We then turned to correlation analyses of 
Hs and water-air temperature difference, Hs and wind speedxwater-air temperature difference, LE 
and wind speed, and LE and wind speedxwater-air vapor pressure deficit for further explanation 
(Fig. 9). Our results suggest that water-air temperature difference explained 45% of the variation 
in Hs. Moreover, the product of water-air temperature difference and wind speed explained 55% 
of the Hs variability. In contrast, about 49% of the variation in LE could be explained by changes 
in wind speed, and 48% of the variation in LE was explained by the product of water-air vapor 
pressure deficit and wind speed. Our results were similar to the findings of Nordbo et al. (2011) 
but opposite with the conclusions of Liu et al. (2009). The reason for the difference is that the 
lake in the study of Liu et al. (2009) is located in a humid coastal area whereas the lake in this 
study is situated in an arid desert. The complex conditions of terrain can produce advection 
(Kaimal and Finnigan, 1994). According to our observations, the atmosphere over the lake 
exhibited a pronounced inversion layer, with lake-breeze circulations. This implies that the 
horizontal transport of latent heat cannot be ignored in arid regions. Meanwhile, salinity could 
affect the heat exchange between lake and atmosphere, while groundwater (temperature of 18°C 
in the study period) recharge supplies a large amount of heat. Therefore, LE and Hs exhibited 
weak run-to-run relationships with the environmental factors (R?<0.50). 
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Fig. 9 Relationships between half-hourly Hs and UxAT (a), Hs and AT (b), LE and U (c), and LE and UxAe (d). 
Hs, sensible heat flux; AT, water-atmosphere temperature difference; U, wind speed; LE, latent heat flux; Ae, 
water-air vapor pressure deficit. 


Several previous studies have shown that the heat fluxes are sensitive to Rn, vapor pressure 
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deficit of atmosphere-lake system and aerodynamic conductance (Williams et al., 1998; Malhi et 
al., 2002; Vourlitis et al., 2002; Kumagai et al., 2005; Loescher et al., 2005; Nordbo et al., 2011; 
Pan et al., 2017), whereas other studies have indicated that the measured LE values are not 
sensitive to these environmental factors (Blanken et al., 2000; Loescher et al., 2005; Liu et al., 
2009). The differences are due to the presence of vegetated surfaces where stomatal conductance 
control is high. 


4.3 Comparison of heat fluxes between lake and land surface 


The variations in Hs and LE observed by the EC system on the nearby land surface (Hu et al., 
2015) were differed from those observed on the lake surface. On the land surface, Hs is strong 
and LE is weak (Hu et al., 2015), whereas on the lake surface, Hs is weak and LE is strong (Fig. 
5). There is a negative correlation between the changes in Hs and LE. In the desert hinterland, LE 
could become negative at times during the year even over the lake, indicating that the latent heat 
is transported downward. This phenomenon has also been observed in other deserts (Wang and 
Mitsuta, 1990, 1992; Harazono et al., 1992). The reason may be due to that the air on the lake 
flows by advection across the dry land and the water vapor flux is in a downward direction on the 
land, implying a downward delivery of water vapor. On the land surface, there were obvious 
changes in Hs during the year (Hu et al., 2015). Specifically, from March to November, Hs was 
mainly positive, meaning that the heat was transferred to the atmosphere. However, from 
November to March, Hs changed dramatically, with a large number of negative values occurring. 
This indicated that the heat was transferred from the atmosphere to the land surface. On the lake 
surface, LE was occasionally negative, implying a downward delivery of water vapor. 


4.4 Evaporation 


The ET values observed on Yindeertu Lake of the Badain Jaran Desert exceeded those observed 
on high-latitude lakes (Schertzer et al., 2000; Rouse et al., 2003; Oswald and Rouse, 2004) but 
were similar to those measured in semi-arid reservoirs (Gallego-Elvira et al., 2010). The average 
ET was 0.6 mm/d during the frozen period for Yindeertu Lake, which is smaller than the 
previously simulated ET (1.3-1.4 mm/d) for Lahontan Lake (Hostetler, 1991). This difference 
might be due to the similar latitude but different longitudes between the two lakes. Furthermore, 
the atmosphere over the lake in our study was unstable at night. As a result, a strong lake-breeze 
circulation formed, which increased the evaporation (accounting for 40% of the total evaporation). 
Because the net solar radiation and water-air vapor pressure deficit are large in the desert and are 
coupled to dry air, evaporation is strong in the desert areas (Zhu et al., 2012). However, ET from 
the lake surface measured in this study is much smaller than those estimated by Chen et al. (2004) 
and Gates et al. (2008a, b). Note that the net evaporation is 10 times larger than the precipitation 
in the study area, indicating that the lake is supplied mainly by groundwater during the study 
period. 


5 Conclusions 


This study measured the radiation components, evaporation rates and heat fluxes by an EC system 
from March 2012 to March 2013 over a lake in the Badain Jaran Desert, China. Both shortwave 
and longwave radiation exhibited seasonal dynamics and large fluctuations. LE exhibited both 
seasonal variations and diurnal dynamics. LE was mainly affected by the changes in wind speed, 
while Hs was primarily affected by the product of water-air temperature difference and wind 
speed. The diurnal LE and Hs were negatively correlated in the open-water period, and the 
effective energy of the lake was mainly directed into LE. Generally speaking, LE was the 
dominant heat flux throughout the year. The average turbulent heat flux over the lake surface was 
greater than the net solar radiation received by the lake, indicating that other energy supplies to 
the lake cannot be ignored. LE showed negative values occasionally both on the land surface and 
lake surface, implying a downward delivery of water vapor. The cumulative evaporation was 
1445 mm/a, and the mean evaporation rate was about 4.0 mm/d. The total evaporation was 10 
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times greater than the total precipitation, indicating that the lake was supplied mainly by 
groundwater during the study period. Continuous and long-term EC measurements will be needed 
to allow an interannual analysis to be performed in the future. 
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